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Quantum oscillations can be used to determine properties of the Fermi surface of metals by 
varying the magnitude and orientation of an external magnetic field. Topological insulator surface 
states are an unusual mix of normal and Dirac fermions. Unlike in graphene and simple metals, 
Berry's geometric phase in topological insulator surface states is not necessarily quantized. We show 
that reliably extracting this geometric phase from the phase offset associated with the quantum 
oscillations is subtle. This is especially so in the presence of a Dirac gap such as that associated 
with the Zeeman splitting or interlayer tunneling. We develop a semi-classical theory for general 
mixed normal-Dirac systems in the presence of a gap, and in doing so clarify the role of topology and 
broken particle-hole symmetry. We propose a systematic procedure of fitting Landau level index 
plots at large filling factors to reliably extract the phase offset associated with Berry's phase. 

PACS numbers: 73.25.-|-i,72.20.Fr, 73.43.Qt 



I. INTRODUCTION 

Topological insulators are bulk insulators with metal- 
lic surface states that are topologically protectedi"— . The 
quantized Hall resistance in two-dimensional topological 
insulators was measured several years ago^. Recently 
measurements have been made of properties of the sur- 
face states of three-dimensional topological insulators. 
Angle Resolved PhotoEmission Spectroscopy (ARPES) 
studies provide direct imaging of near-surface bands^ 
showing dispersions of an odd number of Dirac cones, as 
expected for a topological insulator. Furthermore, these 
band structures are consistent with the presence of a 
non-trivial Berry's phase of tti^ However, these surface 
band structures also exhibit significant band-bending, 
and particle-hole symmetry with respect to the Dirac 
point is broken. 

A complementary method to ARPES to determine the 
Fermi surface properties of a metal is to measure quan- 
tum oscillations such as associated with the Shubnikov 
de Haas or de Haas van Alphen effectsi^ Such experi- 
ments can distinguish between bulk and surface states, 
even if both are gapless (or near-gapless), by tilting the 
magnetic fieldiS. Furthermore, the phase offset (7) of 
the quantum oscillations is related to the Berry's phase 
associated with cyclotron orbits and provides a means to 
experimentally access this important signature of a Dirac 
cone. Experiments on graphen e^^'^^ found evidence of 
the expected non-trivial Berry's phase of tt. There have 
now been a number of quantum oscillation experiments 
on 3D topological insulators J^ii^"— A prevalent observa- 
tion in these studies is that the phase offset is not equal 
to the zero value associated with a Berry's phase of tt. 
In Table U we have listed the different phase offsets ob- 
tained by various groups and materials. The deviation of 
the observed 7 from zero has been variously attributed 
to the Zeeman efFeclji^i^, and the non-ideal Dirac (i.e. 
nonlinear) nature of the surface statesi^ In the present 
work, we reformulate both of these proposals in a precise 



and concrete manner. 
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0.05 ± 0.02 


20 


BiaTeaSe (S4) 


47 


0.32 


20 


BiaTeaSe 


73 


-0.05 


21 


BiaSea 


many 


-0.1 ±0.1 


22 


Bii.996Sno.oo4TeaSe 
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TABLE L The different phase offsets (7) measured for various 
materials and at various chemical potentials (quantified by the 
frequency of the quantum oscillations (Bo), the surface area 
of the cyclotron orbit in teslas). The second last entry was 
performed at more than 20 different chemical potentials. S# 
denotes the sample number in the experiment. 

The aim of this Article is to provide a quantitative 
framework for using quantum oscillation experiments to 
characterise topological insulator surface states. In par- 
ticular, we focus on extracting the phase offset and the 
cyclotron effective mass from experiments. For mixed 
normal-Dirac systems (i.e. Dirac cones with a particle- 
hole symmetry breaking quadratic term), we find that 
there is a key energy scale relevant to quantum oscilla- 
tion experiments: the normal fermion mass multiplied by 
the Dirac Fermi velocity squared, mv^ . The interplay of 
this energy scale with the Zeeman or intrinsic Dirac-gap 
(A) introduces subtleties to quantum oscillations that 
both complicate and enrich the observed phenomena. In 
particular, we find that the phase offset is only quantized 
if the system is either gapless or particle-hole symmetric. 
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In a topological insulator surface state under an applied 
magnetic field, the phase offset is never independent of 
magnetic field because of the Dirac mass gap associated 
with the Zeeman splitting, together with the ubiquitous 
particle-hole symmetry breaking. We outline how one 
can circumvent this difficulty: by fitting the experimen- 
tal data at large filling factors to a simple function (Eq. 
pop ), one can obtain a linearized form which is asymp- 
totically exact in zero field. The intercept of such a linear 
plot as 1/B — 7^ yields the topologically relevant phase 
offset at zero field. We hope that this analysis will be 
valuable to future experimental studies, and will stimu- 
late further studies utilising quantum oscillation experi- 
ments to investigate other topological insulator regimes, 
such as associated with thin films.— 

The structure of the article is as follows. In Sec. II, 
we develop the semiclassical theory of quantum oscilla- 
tions when particle-hole symmetry is broken. In Sec III, 
we introduce a specific model Hamiltonian. In Sec. IV, 
we analyse the quantum oscillations in the model sys- 
tem, and focus on the non-universal phase offsets and 
the subtleties of cyclotron effective mass measurements. 
This section ends with an analysis of expected Landau 
level index plots in these experiments. In Sec. V, we 
propose a robust procedure to extract the phase offset 
from experiment, overcoming the difficulties which led 
to the wide range of values in Table ID In Sec. VI we 
compare our results with existing experiments and show 
that the procedure outlined in the preceding section can 
substantially improve phase offset measurements. In Sec. 
VII, we comment briefly on the relevance of our results 
to spintronic systems. 



II. SEMICLASSICAL THEORY OF QUANTUM 
OSCILLATIONS 



By measuring the longitudinal resistivity of a metal 
sample as a function of external field strength B, one 
observes that it oscillates according to^ 



Apxx oc cos 
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So 
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A. The phase offset in quantum oscillations 

Under a weak external magnetic field, electrons follow 
equi-energy contours about the electronic dispersion. In 
order to construct a semi-classical theory of this situ- 
ation, one usually begins with Onsager's semi-classical 
quantisation condition due to the single-valuedness of 
Bloch wavcfunctions^. The k-space area S{C) of a closed 
electronic orbit C is quantized according to 



S{C)ll=2TT{n + -f{C)) 



(2) 



where is the magnetic length, defined hy Ib = ^ft/ei?, 
and n is an integer. The phase offset 7(C) = 1/2 — 
r(C) /2n, with r(C) the Berry's phase of the orbit (M^, 
given by 



ro(C) 



t J) dk- (uk,a|VkUk,a). 



(3) 



in the band a. This is expressed in terms of the eigen- 
vectors |uk,a) of the fc— dependent Hamiltonian 



H(k) = cxp(— ik • r)iJ cxp(ik • r). 



(4) 



(1) 



where H is, say, a tight-binding Hamiltonian. For nor- 
mal fermions with a quadratic dispersion, or in fact any 
isolated band, T{C) = for any contour C. For mass- 
less Dirac fermions, the Berry's phase is again path in- 
dependent, giving r(C) — TT. Accounting only for the 
bare dispersion contribution to the cyclotron orbit, it is 
tempting to take C to be the equi-energy closed path of 
the zero-field dispersion. Thus 7(C) for a contour which 
includes only the bare dispersion contribution to the en- 
ergy gives a constant 7 = 1/2 for normal fermions, and 
7 = for massless Dirac fermions^. However, for mas- 
sive Dirac fermions the dispersion is no longer linear at 
low energies, and the Berry's phase is no longer path 
independent, (i.e. F = F(C)), and so 7(C) is not, in 
general, quantized<^ 

In a clear and insightful recent article, Fuchs et al^ 
pointed out that for a massive Dirac cone in a magnetic 
field, there is a pseudo-spin orbital magnetic moment 
contribution to the energy that must be taken into ac- 
count for cyclotron orbits {i.e. equi-energy contours). 
The Onsager condition for cyclotron orbits is generalised 
in this case to 



where the oscillation frequency, Bq , is related to the area 
of the Fermi surface enclosed by the cyclotron orbit, and 
7 is a phase offset. This is the Shubnikov-de-Haas effect. 
The de-Haas-van-Alphen effect, which is the oscillation 
of the magnetisation of a sample as a function of B takes 
a similar form, and in two dimensions is given by Eq. 
p4)) . The phase offset, 7 is related to the Berry's phase 
of the cyclotron orbit 1^ A key issue we clarify in the 
present work is the role of the Berry's phase in the phase 
offset 7. 



5(e)Z| = 27r(n + 7„(e)). (5) 

The quantity jai^) is not, in general, equal to 7(C), since 
it also includes the pseudo-spin magnetisation compo- 
nent. It is 7a (e) that gives the relevant phase offset, and 
not 7(C) in Eq. ©. 

Somewhat surprisingly, for the particle-hole symmetric 
massive Dirac cone, this extra contribution to the phase 
offset due to the pseudo-spin orbital magnetic moment 
exactly cancels the energy-dependent part of the Berry's 
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phase, such that 7a (e) = after all<^ By inverting Eq. 
(O, it can be shown that the semi-classical Landau levels 
exactly reproduce the fully quantum mechanical ones in 
the specific case of a massive Dirac cone. For this reason, 
Fuchs et al. label 7a (e) (7l in their nomenclature) the 
'Landau index shift', and they associate it with a winding 
number of the orbit - a topological quantity. Therefore, 
the phase offset is strictly 7=1/2 for normal fcrmions, 
and 7 = for particle-hole symmetric Dirac fcrmions. 



B. A general expression for the phase offset in 
quantum oscillations: broken particle-hole symmetry 



A further goal of this work is to generalise the contri- 
bution of Fuchs et al. to systems with broken particle- 
hole symmetry. This extension applies to many systems, 
most notably at present, topological insulator surface 
We find that in the case of broken particle- 



.34.35 



statei 

hole symmetry, the winding number can no longer be 
associated with the magnetic moment-adjusted Berry's 
phase, and so the strict quantisation of 7Q(e) is removed. 

Consider a general two band system. Following Fuchs 
et al.^ the energy of an electron in band a is 



«5(e„(k))/| + i(ea(k)-ea(k)) 



dVc,{e) (10) 
*Q(k) ' 



In the last equality we made use of Eq. ^ and Eq. ([9|) . 
Rearranging the above expression to obtain 5'(eQ,(k)), 
and recognising that we already know S'(eo,a(k)) from 
Onsager's condition, Eq. ([5]) we obtain 



7a (e) 



1 

2^ 



(11) 



This result for the phase offset in particle hole symme- 
try broken mixed Dirac-normal fermion systems directly 
leads to Eq. the central result of this paper. Two 

limiting cases are worth mentioning. Firstly, for particle- 
hole symmetric Hamiltonians, eg, a = — eo,a, and the sec- 
ond and third terms on the right together reduce to 
a multiple of the winding number already obtained by 
Fuchs et al., which is simply 1/2, giving ^a{^) ~ 0. Sec- 
ondly, for a degenerate two band system, eg, a — ^o.a, and 
r = 0, trivially giving 7a(e) = 1/2. In general, however, 
no such simplification is possible. 



e„(k) =£o,a(k)-X„(k)-S, 



(6) 



where eo,Q(k) is the bare dispersion energy in zero field, B 
is the external magnetic field, and the pseudo-spin orbital 
magnetic momen t^^'^^ 



A^„(k) = -((f-rc)xj) 

(7) 

= |-(eo,a(k)-eo,a(k))l)o(k), 

where Vc is the centre of mass position of the wavepacket, 
j is the current operator, and a denotes the band which 
is not a. We emphasise that Eq. ([7]) is a property of 
Bloch electrons in a magnetic field, and has nothing to 
do with spin, fla (k) is the Berry curvature which is the 
local quantity given by 



Jla(k) = Vk X i(uk,a|VkMk,Q) 



(8) 



and is related to the Berry's phase by Stokes' theorem 
such that 



r„(c) = / r!„(k)dk, 

Js 



(9) 



where S is the area enclosed by C. Eq. ([7]) generalizes 
Appendix C in Ref.([29f) to the case of broken particle- 
hole symmetry. Including the magnetic moment in Eq. 
([5]). wc can obtain the full quantisation condition for a 
generic two-band model, by finding 



C. Semi-classical theory of the 
de-Haas— van-Alphen effect 

The Lifshitz-Koscvich theory of the de-Haas van- 
Alphen effect for a single species of free electrons has 
been obtained by Champel and Minee^*^. Starting from 
an expression for the Green's function of an electron in 
a magnetic field, they obtain the density of states, and 
then the thermodynamic potential. From here one can 
directly obtain the magnetization, and see that it oscil- 
lates as a function of magnetic field. In particular, this 
demonstrates the direct relation between oscillations in 
density of states as a function of field strength, and how 
they directly lead to oscillations in magnetization. 

By inverting Onsager's condition, Eq. (O, in terms 
of energy levels, and thus replacing the energy levels 
in a magnetic field with an expression for the quan- 
tized orbits 5(e) and the phase offset 7, Luk'yanchuk 
and Kopelevich^ generalised the treatment of the de- 
Haas- van-Alphen (dHvA) semiclassically, for arbitrary 
systems. 

In a magnetic field, the electrons in band a undergo 
cyclotron orbits with frequency 



Wc(ec<) 



eB 



m*(e)' 

where the cyclotron effective mass is 

1 dS^e) 



m;(e) 



27r dca 



(12) 



(13) 
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The 2D result obtained by Luk'yanchuk and 
Kopelevich^^ of which we are interested here, is given 

by 



'^dSaie)/de ^ sinh(AO^ 



X sin 2Trl 



2TTheB 



7a (e) 



cos 27r/ 



gs^J■BB 



where 



A 



2n^kBTdS{e) 



de 



(14) 



(15) 



and h/r is the Landau level broadening. 

Quantum oscillations provide an experimental method 
to determine 7a (e). By assigning integers (n) to the oscil- 
lation maxima/minima in the magnetisation or resistance 
as a function of inverse field, one can extrapolate these 
to 1/_B — ^ to obtain 7c((e) which is the n— intercept. 
Modulo 1. 



III. TOPOLOGICAL INSULATOR SURFACE 
STATES 

A. Model Hamiltonian 

Consider a two-band Hamiltonian, where the basis is 
not necessarily spin, but can be thought of as a pseudo- 
spin, and may be spin, sub-lattice, on-site orbitals, etc, 
and which contains both normal fermions, and massive 
Dirac fermions, given b y^s, 34^35, 37 



H 



2m 




hvp{ky + ikx) \ 
ik^) -A ) 

(16) 

where there are four free parameters, m, /z, VF^ and A. 
m is the mass associated with the normal part of the 
spectrum, ^ is the chemical potential, vp is the velocity 
of the Dirac part of the spectrum, and A is the Dirac 
gap. 

This Hamiltonian, at A = 0, is the low energy surface 
theory of a time reversal invariant topological insulator. 
Following the classification scheme outlined by Schnyder 
et al^, we see that the quadratic term breaks particle- 
hole and chiral, or sub-lattice symmetries. Since we are 
concerned with electrons then, time reversal is antilin- 
ear, and the subclass of the topological insulator is AH, 
and the topological number is a Z2 invariant. Generi- 
cally, finite A gives a Dirac mass gap, and so the three 
dimensional system is no longer a topological insulator. 

In the absence of an orbital magnetic field there are 
two bands, with dispersion 



«(k) 



2 7,2 



2m 



A2 



(17) 



where a = ±. 

The Landau level spectrum can be calculated exactly, 
givmgii* 



= hojQn ± 



' 2hvj,eBn 



hujQ 



gsf^BB \ 
2 J ' 



(18) 



where luq = eB/m. In the limit vp — > 0, we obtain the 
usual normal fermion Landau level spectrum, and in the 
limit TO — )■ cxD, we obtain the Dirac fermion spectrum: 



^ I f} Wo{n + ^) 
I y/2hvj,eBn 



(Normal fermion) 
(Dirac fermion) 



(19) 



which arc quoted at A = 0. Respectively, the phase offset 
for quantum oscillations can be 'read off', as 7 = 1/2 for 
normal fermions, and 7 = for Dirac fermions. In the 
mixed case of Eq. p^ . however, there is no constant 
which can be so readily extracted from the Landau level 
expression. 

We can define a useful 'interpolation parameter', given 
by 



mvt 



(20) 



which for fj,,m > gives < 77 < 1. The interpolating 
parameter gives as two limiting cases at finite fi: 




(Normal fermion) 
00 (Dirac fermion). 



(21) 



Generally, rj small gives a Rashba spin orbit coupled 
2DEG)^ where vp is the Rashba spin splitting, and A = 
gs^^sB is the Zeeman splitting ; 77 « 1 gives, for example, 
graphene at A = 0, and boron nitride at A 7^ 0,— 1^ 
and < 77 < 1 gives the surface states of topological 
insulators^^^, where A can be the Zeeman splitting or 
a tunnel splitting for thin films^. Although we are most 
interested in the latter regime, our results are generally 
valid for all values of 77 and A. 



B. Broken particle- hole symmetry 

In real materials topological insulator (TI) surface 
states are not simply described by massless Dirac cones. 
ARPES experiments show surface states with significant 
band-bending and broken particle-hole symmetry with 
respect to the band crossing pointi^i^i^. For this reason, 
Hamiltonian Eq. is the relevant low energy Hamilto- 
nian to describe topological insulator surface states. This 
has been confirmed by symmetry argumcnts^^^. In Ta- 
ble [ni we show some of the values for the mass to, and 
the Fermi velocity vp characterising the Hamiltonian Eq. 
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FIG. 1: Model band structure with particle-hole asymme- 
try, highlighting the effect of the /2m term relative to 
the chemical potential. In a particle-hole symmetric system, 
e+(fc) = — e_(fc). At finite m, away from fc = 0, the chemical 
potential is shifted up for either band with respect to the sym- 
metric case, by an amount S^/n ~ ij,/{2mvp) = (l/?7 — l)/2. 
In the limit /i — or m — >■ oo, the shift 5jj,/jj, — >■ 0, restoring 
particle-hole symmetry. 




FIG. 2: For all energies, there are two distinct cyclotron orbits 
(Fermi surfaces). However, for real materials the Hamiltonian 
Eq. p6|) is an effective one for low energies, and the outer 
orbit may not actually exist. Nevertheless, for systems with 
small ?7, such as Rashba systems,— there are two possible 
orbits, indexed by /3 = ±1 for outer and inner orbits. In the 
figure, ^ > 0, so the inner orbit (/3 — —1) circulates around 
the a — 1 band, and the outer (/3 — 1) orbit circulates around 
the a = — 1 band. For fi < 0, both orbits circulate around 
the a = — 1 band. 

< rj < I, the system is neither purely Dirac nor purely 
normal, but is a mixture of the two. In particular, near 
k = 0, for any finite rj, the electrons are Dirac fermion- 
like, and far away from fc = 0, they become normal 
fcrmion-like. We can rearrange Eq. (|2ip to quantify this 
statement in terms of fi 



measured by ARPES. It is clear from this data that 
for typical doping levels, the topological insulator surface 
states are truly mixed normal-Dirac systems, and cannot 
be considered to be approximating either regime. 



Material 


m/nie 


vf (ms 


jj. (meV) 


V 


Ref. 


Bi2Se3 


0.25 


5.0 X 10^ 


300 


0.54 


42 


BizTeaSe 


0.13 


3.4 X 10^ 


130 


0.39 


i5 


BiaTes 


3.8 


3.9 X 10^ 


250 


0.93 


43 


Bii.sSbo.sTei.ySei.a 


0.32 


4.6 X 10^ 


140 


0.73 


16 



TABLE II: Estimates of model Hamiltonian parameters for 
several 3D topological insulator materials. Normal fermion 
mass m and Dirac fermion Fermi velocities vp are estimated 
from ARPES. At the values of the chemical potential jj, quoted 
in these works, we have also calculated the parameter rj, de- 
fined in equation Eq. (|20|) . 

The effect of broken particle hole symmetry on the 
transport properties of topological insulators has been 
investigated theoretically;^ and experimentally!^ 

We can think of the two limiting cases of Eq. (PT|) 
in terms of particle-hole symmetry. In particular, pure 
Dirac fermions have particle-hole symmetry for all fi, 
whereas pure normal fermions have two degenerate bands 
and so have broken particle-hole symmetry for all fi. For 



We can say that r] provides a measure of broken particle- 
hole symmetry with respect to fc = 0. The effect of bro- 
ken particle-hole symmetry on the band-structures is elu- 
cidated in Fig. [H 

C. Gap opening: Zeeman and tunnel splitting, and 
beyond 

Mikitik and Sharlai have argued that, in general, for 
Bismuth-based 3D topological insulators, 7 = in the 
weak field limit 4^ It has been shown that the symmetries 
of the Bismuth based TIs don't allow for a non-zero gap A 
in Eq. p^| Mi^. and so apart from the Zeeman splitting, 
and as we shall show, 7 is strictly zero. 

However, it is clear from Table U that 7 is rarely ob- 
served to be zero, but in fact deviates substantially from 
zero. This implies that the g-factor may be 10 — 100 
in topological insulator surface statesi^^. This cor- 
responds to a Zeeman gap of 6 — 60 meV in a lOT 
field, which, when compared to typical Fermi energies 
of 50 — 300 meV (Table |l|, is not insignificant. 

In topological insulator thin films, the surface states on 
the two faces of the sample hybridise, and thus open a 
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FIG. 3: (Color online) The phase offset 7+ as a function 
of A, the gap parameter, and rj, the dimensionless param- 
eter tliat tunes between normal fermions (77 = 0) and Dirac 
fermions (77 = 1). It is this quantity, and not simply the 
Berry's phase, that determines the phase offset in quantum 
oscillation experiments. If either particle-hole symmetry is 
maintained {rj = 0, 1), or the system is gapless, the phase off- 
set is a constant. If neither condition applies, the phase offset 
is a continuous function of 77 and A. The solid lines denote 
constant values of 7+. 



mined from the full dispersion Eq. p7|) . and the Berry's 
phase at fixed energy is implicitly dependent upon the 
full Hamiltonian. Moreover, the a = — 1 band has both 
an inner and an outer orbit, as shown in Fig. [5J In par- 
ticular, solving the dispersion at fixed energy e for k^, we 
obtain 

kj = |j(e + mvl + (3^A^ + 2emvl + m^vj,), (24) 

where /3 ~ ±1. The Berry's phase and phase offset 7 
will then acquire an inner/outer orbit index. For a = 1, 
/3 = — 1 always, as there is no outer orbit in this case. 
However for a ~ —1, and e < 0, (3 ~ [—1, 1] coexist, and 
for e > 0, /3 = 1 always. At fixed energy then. 



(25) 

With the Berry's phase known, we are able to find the 
phase offset for quantum oscillations from Eq. (|TT|) . For 
the inner orbits (/3 = —1), the relevant case for topolog- 
ical insulators, we obtain 



So even in Bismuth-based TIs, a non-zero A in Eq. 
([in]) can occur. In fact, thin films below the critical thick- 
ness (6 quintuple layers^) have already been successfully 
produced in the lab, and the gap observedi^i^— . 

There are other mechanisms for gap formation which 
have been suggested, and even realised, including exci- 
tonic pairing gaps^ii^ and a possible Higgs-type mech- 
anism near the phase transition from the topological 
insulating phase to the topologically trivial insulating 
phasci^ The formation of a gap is also crucial to var- 
ious topological effects becoming manifest, including the 
as-yet unverified topological magneto-electric effect^. 



7a.- (e) 



aA 



1 



_2c 

mv'i 



(26) 



This is the central result of this paper, which shows 
that the phase offset expected in quantum oscillation ex- 
periments on topological insulator surface states is non- 
universal, but in fact depends upon the energy scales of 
the gap, and a mixed normal-Dirac fermion energy scale, 
mvp. 

For the outer orbits, only relevant for the a = — 1 band, 
and, for example, Rashba systems where is small^^ 



IV. QUANTUM OSCILLATIONS IN 
TOPOLOGICAL INSULATOR SURFACE STATES 

A. Non-universal phase offset in quantum 
oscillations 

The Berry's phase is given by Eq. ([31), and is readily 
calculated for the Hamiltonian Eq. ([T5|. In particu- 
lar, noting that the Berry's phase is a function of the 
wavefunctions only, the k^ term in Eq. does not 

contribute for a path at fixed momentum space radius k, 
and so at all 77, we have 



7-,+ (e) = 7- -(e) + 



(27) 



Where x = -^/A^ + 2emvp + m?Vp. We have included 
this result for completeness, but from now on will focus 
on the inner orbits, relevant for topological insulator sur- 



face states, and drop the /3 index in Eq. 

There are three particularly interesting limiting cases: 



7a (e) = <^ 0, 
0, 



vp (Normal fermions) 
771 — > cx) (Dirac fermions) 
A ^ (Mixed gapless) 



(28) 



^a{k) = 7ra(l 



v/A2 + h^i 



(23) 



which is the same result obtained in Ref.[29| for the 
case 7Ti — > 00. At fixed energy however, fc(e) is deter- 



The second and third limit together are particularly in- 
teresting, as they imply that for infinitesimally small vp, 
at A = 0, 7 = 0, and so 7 jumps discontinuously from 
to 1/2 at vp = 0. This discontinuity is apparent in the 
vertical line of Fig. [3]at A = as r/ — > 0. 
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The second limit gives the same result obtained previ- 
ously by Fuchs et al.^: in the presence of particle- hole 
symmetry and a Dirac gap the phase offset is zero, even 
though the Berry's phase F is not equal to tt and varies 
with the chemical potential. 

Fig. [3] shows the expected phase offset for the elec- 
tron band of gapped systems with both normal and Dirac 
components in terms of the two parameters A = A//i and 
T]. The phase offset Eq. ([^5]) can be written in terms of 
these as follows: 



7„(A//i,7;) 



a(l-77)A^ 



2v/r;(2-r;) + (A/^)2(l-r;)2 



(29) 



The phase offset varies as a function of the energy scale 
mvp. The horizontal line t] ~ corresponds to purely 
normal fcrmions, and 77 = 1 to purely Dirac fermions. We 
can see that along these two lines, 7 is indeed quantized 
to 1/2 and respectively, as expected. The topological 
result of Fuchs et al^ can be seen along the horizontal 
77 = 1 line of the upper panel. Clearly in the more general 
case, the same topological arguments made in that case, 
relating 7 to a winding number, no longer apply, as 7 is 
no longer quantized. Moving away from the horizontal 
line at 77 = 1 corresponds to the breaking of particle- 
hole symmetry. Moving away from the line A = can 
correspond to breaking inversion or time reversal sym- 
metry, depending on the basis, if the low energy theory 
Eq. (jl6p is applicable at a high symmetry point in the 
Brillouin zone. In this particular regime, and considering 
just a single band, we recover the result of Mikitik and 
Sharlai,— that 7 = irrespective of the extent of particle- 
hole asymmetry. However, in more general systems, this 
does not correspond necessarily to the breaking of any 
particular symmetry. A prominent example is graphcne 
with a staggered sub lattice potential, where the symme- 
try being broken is a sublattice one^. 

So we see that in the general case, if either of the two 
symmetries are retained, such that A = or 77 = or 1, 
7 is a constant, and independent of energy. If both are 
simultaneously broken, then the phase offset is a function 
of the gap (external field), the chemical potential, and the 
normal- fcrmion-Dirac fermion energy scale mvp. 



,+mv'^p{ 1+13 



A2 



2e. 



+ 1 



(30) 

as outlined in Fig. [2j 

We can easily check that as — >■ 0, we obtain the 
normal fermion result of S"(e) = 27r77T,, and by expanding 
to all non-zero orders as 777, 00, we obtain the usual 
Dirac fermion result of <S"(e) = ire/vp. 

However, we note that in the mixed case, experimen- 
tally determining Ec^. pop docs not directly give m ot vp, 
but contains an expression with a pair of energy scales 
mvp, and A. In fact, for a mixed system such as this, we 
obtain from Eq. p^ . X (x ^ = 2'!Tm* , and so we have 



= 1 - 



VA^ 



2emVp 



(31) 



In Fig. |4]wc have shown the ratio of the cyclotron mass 
to the normal fermion mass, the relative cyclotron mass, 
given by 



771 

m 



A2 



(1 - r;)2 + 77(2 - 7j) 



(32) 



This figure shows the gap and chemical potential depen- 
dence of the cyclotron mass, highlighting that decreasing 
the gap or increasing the chemical potential tunes the cy- 
clotron orbit to a more Dirac-like part of the spectrum. It 
also demonstrates that the interpolation between normal 
and Dirac fermions smoothly interpolates the relative cy- 
clotron mass between 1 and 0. Note that the vanishing of 
the relative cyclotron mass at 77 = 1 is due to 771 ^ 00. In 
the Dirac limit, the cyclotron mass is 777,* = ^/vp. Nev- 
ertheless, at either 7/ = or 77 = 1, the effective mass 777* 
is gap and chemical potential independent. In between 
however, this is not the case, and the effective mass is a 
function of the gap and the chemical potential. 



C. Quantum oscillations and Landau level index 
plots 



B. Cyclotron effective mass 

Using the temperature dependence of the magnitude 
of the magnetisation at fixed field, one can fit the pref- 
actor in Eq. that is A/ sinh(AZ), where A is given by 
Eq. . Usually, one can then determine the quantities 
777 (normal fermions) or vp (Dirac fermions). However, 
in our mixed system, the cyclotron effective mass is more 
complicated, being given by Eq. (jl3|) . where, reintroduc- 
ing the inner-outer orbit index. 



The oscillations observed in resistivity (Eq. ([T|)) and 
magnetisation (Eq. (|14p ) experiments have extrcma at 
completely filled or empty Landau levels, to which we 
can assign integers and half integers 77 (77 1/2). There- 
fore a plot of the location of the minima/maxima as a 
function of magnetic field can be used to identify the 
magnetic field values at which one has a filled/ half-filled 
Landau level. From Eq. ([1]) and Eq. (|14p . it is clear a 
sensible ordinate for a plot of minima/maxima is the in- 
verse field strength, 1/ B. To a good approximation then, 
{i.e. only taking the first term in the scries in Eq. (|14p ). 
the condition for integer filled Landau levels is 
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FIG. 4: (Color online) The relative cyclotron mass vn* jm as 
a function oi A./^ and band bending. At 77 = 0, m* = m. At 
?7 = 1, m — >■ 00, so m* /m — >■ 0. As the gap is decreased, or the 
chemical potential increased, the system becomes more Dirac- 
like. At 77 = or 77 = 1, the cyclotron mass is independent of 
the gap/chemical potential, however in between the two this 
is not the case. 



2TTheB 



7a (e) = n - A, 



(33) 



where A ~ 1/2(0) for the minima (maxima) of the lon- 
gitudinal resistivity, and A = 3/4(1/4) for the minima 
(maxima) of the magnetization. Further, we note that as 
1/ B — > 0, the intercept of the index plot with the n— axis 
(offset by A) gives 7(e). 

In order to determine the index plot, a knowledge of 
both the area of the orbit, S'(e), and the phase offset 7(e) 
is required. We have already given 7(e) in Eq. (j26p . and 
plotted 7 for all 77 and a range of gap parameters in Fig. 
[3l S{e) is given in Eq. pop . however it is interesting to 
note the limiting cases: 



S'a(e) 



2iTme/h^ , vp ^ (Normal fcrmions) 
^{6^ ~ A'^)/h'^Vp, m 00 (Dirac fcrmions) 



(34) 

Firstly, we note that 7, Eq. (pS)) . only varies as a func- 
tion of the gap parameter in mixed normal-Dirac systems, 
being constant with respect to the gap in purely normal 
or purely Dirac systems. Therefore, for purely normal 
or purely Dirac systems, one expects the extrapolation 
of an index plot as 1/B to give 7(e) = or 1/2. 
Second, wc notice that for gapless Dirac fermion and all 
normal fermion systems. S{e) does not vary with the ex- 
ternal field; and so the index plot is strictly linear. For a 
Zeeman gapped massive Dirac fermion system however, 
there is a nonlinearity in the index plot due to the Zee- 
man term A = gs^sB. Fortunately for a material such 
as graphene, there is no Zeeman gap, and so the index 
plot reliably gives the expected intercept for relativistic 
fcrmionSfii In all these cases then, correctly extrapolat- 



< 
I 



* Extrema 
Fit 

Asymptotic limit 



0.05 



0.1 



0.15 



0.2 



1/B (T"') 



FIG. 5: (Color online) Robust determination of 7b^o from 
an index plot. The crosses represent the extremum of the 
resistivity or magnetisation from Eq. ^ or Eq. p4p corre- 
sponding to filled Landau levels. The solid line is the small 
field fit, Eq. (|36p . It is clear that even for a system with 
pronounced Zeeman splitting (here gs = 50), the fit is reli- 
able at all filling factors. The dashed line is the linearized, 
asymptotic fit {d^/dB = in Eq. ((36])), which is asymp- 
totically exact as n — > 00. Extrapolating this line back to 
l/B — >■ gives the topologically relevant phase offset, 7s^o, 
which shows the Dirac nature of the surface states. System 
parameters are relevant to Bi2Te2Sei^: vf = 3.4 x 10^ms~^, 
m = 0.13me, Qs = 50. 



ing the intercept at 1/B yields a constant field and 
energy independent value of 7. 

In a system with non-zero Zeeman splitting, and both 
normal fermion and Dirac fermion terms however, 7(e) 
becomes field and energy dependent, as does S{e). When 
calculating a Landau level index plot from the max- 
ima/minima of the oscillations (Eq. in mixed sys- 
tems such as topological insulator surface states then, the 
index plot as a function of becomes uniquely nonlin- 
ear according to the energy scales of the Zeeman splitting 



and, as is clear from Eq. ([3D]) and Eq. (|26p . 



In Fig. 



Othis nonlinearity with increasing field can be clearly ob- 
served, and has been pointed out elsewhere^^^. 

It is clear from the above discussion that extrapolat- 
ing an index plot to \/B yields an estimate of 7(e). 
However, from Eq. ([26]) . we see that 7(e) in the case 
of mixed normal fermion-Dirac fermion systems depends 
on the energy scales A and mvp. For this reason, the 
index plot cannot possibly discriminate between m and 
vp. The cyclotron effective mass (Eq. (|13l) ). on the other 
hand, intrinsically bears such a distinction. Through a 
combined approach of measuring the temperature depen- 
dence of the oscillation amplitudes, and thus obtaining 
TO* (Eq. ([T3)) ). together with index plot measurements 
to determine the energy scales mvp and A, one can de- 
termine both the normal fermion mass m, and the Dirac 
fermion velocity vp. 
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USING LANDAU LEVEL INDEX PLOTS TO 
EXTRACT THE PHASE OFFSET 



Taking into account the effect of a bulk Fermi 
surface 



The quantum oscillation experiments performed thus 
far— li^"— , tend to report a substantial bulk contribution 
to the measured resistivity. Although there is a bandgap 
at the r point of the surface Brillouin zone, as clearly 
observed with ARPES^ii^, the bulk 3D Brilluoin zone is 
not usually gapped throughout, but has a Fermi surface. 

As our Hamiltonian, Eq. ([T5)) . is merely the effective 
surface theory of a bulk system, we must take into ac- 
count the effects of the bulk on quantum oscillation ex- 
periments. In particular, we must consider the following 
expression; 



Ne = 



Pb{e)de 



Ps{e)de 



(35) 



where iVg is the number of electrons, which is a con- 
stant, pb is the density of states of the bulk bands, and 
Ps is the density of states of the surface bands. As we 
switch on a magnetic field, remains constant, Ps(e) 
thus changes due to the opening of a Dirac gap, and thus 
p varies. Since the surface band is a two-dimensional 
band, of which there is only one, whereas the bulk bands 
are three-dimensional and there are many, the number of 
bulk carriers is much greater than the number of surface 
carriers. 

Therefore, there are two distinct regimes. 

1. If there is a bulk Fermi surface, the contribution to 
Ne will come almost entirely from the first term in 
Eq. ([55)) . and so the chemical potential of the sur- 
face band will be dictated by the change in chemical 
potential of the bulk band. In the limit where the 
bulk Fermi surface shifts by a negligible amount, 
this is equivalent to taking the surface theory, Eq. 
(|16p . in the grand canonical ensemble, where p re- 
mains constant, and the Fermi wavevcctor kp, and 
therefore S{p) varies. 

2. // the hulk is completely gapped, the first, bulk, 
contribution to Eq. p5p remains constant, and 
the change in p is dictated entirely by the sur- 
face theory. This corresponds to the case most 
frequently considered, where it can be shown that 
the Fermi wavevector remains constant (Luttinger's 
theorem^), and thus S{p) remains constant as p 
varies. 

Xiong et. al^ have pointed out a second subtlety in 
the presence of a bulk Fermi surface. They argue that 
when measuring the longitudinal resistivity, the contri- 
bution of the bulk resistivity demands that filled Landau 
levels be associated with resistivity maxima, rather than 
resistivity minima. 



B. Extracting the phase offset at small fields 

Experimental determination of the index plot in this 
case of mixed Dirac and normal fermions is, as we have 
discussed, fraught. In the presence of an unknown gap, 
an intercept with n yields a non-quantized value of 7 
which can sit anywhere on an equal-7 contour in Fig. [31 
Furthermore, if the gap is from Zeeman splitting, it is lin- 
early dependent upon B, and so, as pointed out by Taskin 
and Ando2^ as well as Seradjeh, Wu and Phillips^^, and 
shown in Fig. [51 will not produce a straight line in the 
index plot as a function of inverse field strength. 

As the phase offset 7(e) is only non-constant when the 
system is both gapped and the electrons arc mixed nor- 
mal and Dirac fermions, we can use this fact to determine 
whether the system has a gap at zero field. It is interest- 
ing to note that due to both of these stipulations being 
simultaneously required, this ability of quantum oscilla- 
tions to extract an intrinsic gap is only possible in the 
unique mixed normal-Dirac systems we are considering 
here. 

In the case of a system with no intrinsic band-gap, we 
would like to measure jb~>o = 0, the 'topological' result 
for Dirac fermions. If the system has a tunnel-split gap, 
or some other more exotic gap, then "fB=o 7^ 0- 

In order to circumvent the difficulties outlined above, 
and in order to systematically obtain a reliable value for 
the topological phase offset, or a measure of the intrinsic 
band-gap. we propose the following procedure. Starting 
with Eq. (|33p . we make a small field expansion and re- 
arrange to obtain 



A 



Bo „d7 
-B-^^^'-^dB 



B, 



(36) 



B->0 



where 



Bn = 



2'nhe 



(37) 



is the area enclosed by the cyclotron orbit with no Zee- 
man correction, and 



d7a(e) 



dB 



agsPB 



2e 
mv'i 



1 



(38) 



The constant C is, in general, unknown. In case 1 out- 
lined in the previous subsection, where there is a bulk 
Fermi surface and the change in p with increasing mag- 
netic field is negligibly small, S{e) varies as a function 
of magnetic field. The variation of 5(e) with magnetic 
field can be approximated by expanding Eq. pop with re- 
spect to the field. The second term in such an expansion, 
which goes as B^, is oc dj /dB\B-^Q- In this case then, 
C has a contribution from the gradient of 7, as well as 
a contribution from the change in the area of the Fermi 
surface at p, which is proportional to the gradient of 7 
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and the ratio of the Zeeman energy (A^ = gsfisB) to the 
cyclotron energy of the normal fermion (ftwo = heB/m). 
In case 2, S{e) is a constant, and so C = 1. Therefore 




1 + -^^j /i constant (bulk Fermi surface), 
1, S{e) constant (bulk gap). 



(39) 

The topologically relevant quantity is 7B_i.o, which is the 
intercept of Eq. (|36p with the n — A axis in the limit 
d-y/dB 0. 

Our proposal is thus as follows: 

1. Plot the experimentally obtained extrema in the re- 
sistivity or the magnetisation as a function of 1/B. 

2. Fit the data to the fitting function 



n-A=-^ + Ai+ A2B, (40) 

where Bq, Ai, and A2 are constants. 

3. The asymptotic low field limit {B — > 0) is equiv- 
alent to ^2 ^ and yields a straight line, whose 
intercept with the n axis is 75^0- Therefore, the 
topologically relevant phase offset is Ai. If = 0, 
the surface states are gapless, and contain a Dirac 
component with a Berry phase of tt. If Ai ^ 0, 
then the system may have an intrinsic band-gap. 

Eq. pop , together with theoretically obtained minima 
of Eq. dl]) or Eq. (HH), is plotted in Fig. [5] for typical 
parameter values. The fitting function is seen to be a 
good fit for all n. The topologically relevant asymptotic 
form of Eq. has also been plotted in Fig. [SJ where 
7B-i.o — 0, as expected for topological insulator surface 
states with no intrinsic band-gap. 

We emphasise that producing a straight line fit to the 
experimental points, rather than following the procedure 
outlined above, will only yield a reasonable estimate of 
7B_i.o for very large n. As can be seen from Fig. [Sj the 
asymptotic approach of the data to the straight line fit 
is slow (approaching as l/i?~^). A more reliable method 
is to fit Eq. (|40l) to the data, and in this way extract the 
topologically relevant zero field phase offset. 

VI. COMPARISON WITH EXISTING 
EXPERIMENTS 

As has been pointed out by Taskin and Ando^^, the 
nonlinearity of the index plots for large B is evident in 
many existing studies, and the extrapolation to l/i? — > 
by fitting a straight line through the data points consis- 
tently yields 7 ^ 0, 1/2 (see Table U). It is also possible 
that for certain values of ^s, or external fields, the index 
plot will curve so much that the extrapolation spuriously 
yields 7 = 1/2 for a topological insulator surface state. 




1/B (T-^) 

FIG. 6: Determination of the zero field phase offset 7b_»o 
using the fitting function Eq. (|40|l (solid lines) on experi- 
mental data for Bi2Te2Sei^, and Bi2Sejr^. For Bi2Te2Se, we 
liave shifted n index by 1/2, following Xiong et ali^, and 
obtain 7s^o ~ —0.1, whereas a linear fit to the data yields 
7s->o = —0.78. For Bi2Se3, we obtain 7s->o = 0.1, whereas 
a linear fit yields 7b-»o = —0.36. The dashed lines are the 
asymptotic plots of Eq. (|36p . with A2 — 0, whose n — A 
intercept gives the topologically important zero field phase 
offset. 

In fact Analytis et ali^ obtained such a result with their 
sample 2 (see Table [l|. It is clear from our results that 
these findings are to be expected, and the problem may 
be systematically circumvented by fitting Eq. (|40l) to the 
low field data and and thus obtaining a reliable estimate 
of 7b->0i outlined above. 

In Fig. [5] we have performed just such an analysis 
for two existing experimental studies on Bi2Te2Se>i^, and 
Bi2Se3i^. Following the discussion of Xiong et al^, we 
have shifted the index identification n in the former by 
1/2. The reason for this being that the bulk contribution 
to the longitudinal resistivity causes the maxima, rather 
than the minima, to be associated with a completely filled 
Landau level. If this is the case, then a linear extrapo- 
lation of the data yields an intercept of Jb^o = -0.78. 
On the other hand, using Eq. (PU)) . our fit becomes 



n--w h0.1-0.05B (41) 

2 B 

Comparing with Eq. ([55)) then, we find Bq = 60 T, 
in agreement with the value obtained in the original 
paper—, and 73^0 = 0.1, which is close to the expected 
value of 0. From Eq. (|55)) . and the expression for C at 
fixed fjL, Eq. p9p . we see that if we use known values of m 
and vp, we can estimate the g- factor from the third term 
in the above result. Using this value we obtain gs ~ 60 
for 77 = 0.39, which is in broad agreement with similar 
estimates given elsewherei^i2£. 

Turning to the case of Bi2Se3i^, we obtain 



11 



28 

'b 



-0.1 -0.018S 



(42) 



Again, we obtain a good agreement with Bq = 28 T, and 
a zero phase offset of ^b^o = 0.1, whereas a hnear fit 
yields 7b-s.o = —0.36. Using the material parameters of 
Bi2Se3 from Table [TTl we obtain broad agreement with 
estimates given elsewhere, namely gs ~ 50 for -q ~ 0.54. 

This brief comparison with experiment highlights the 
need for further experiments on higher purity samples 
with smaller Dingle temperatures. This will enable ex- 
tending the index plots to lower fields. 



VII. APPLICATION TO SPINTRONICS 

For 77 small, the Hamiltonian Eq. corresponds to 
a typical spintronics system: a two dimensional electron 
gas, such as GaAs and InAs quantum wells with Rashba 
spin-orbit interaction)^ For vp small, the two bands be- 
come nearly degenerate, and both contribute to quantum 
oscillations, i.e. the sum over a must be completed in Eq. 

A robust treatment of these systems is not within the 
scope of the current work, as inter-band effects contribute 



a significant third portion to the oscillation 



,54.55 



on top 



of the oscillations of the two intra-bands, and so a third 
term is required in Eq. ()14|) . as inter-band tunnelling 
becomes possible. 

However, the results of the present work do allow one 
particularly striking prediction for spintronic systems. 
According to Eq. ([26]) . and Fig. [3l for any finite Dirac 
velocity, or in this particular case, Rashba spin-orbit in- 
teraction, as A 0, 7 0. For very small Rashba terms 
though, this transition from 7 « ±1/2 to 7 = becomes 
almost step-like. Therefore, for all but the tiniest of Zee- 
man fields, one might expect to measure 7 = 1/2, even 
though the B ^ limit indeed gives 7 = 0. For low 
fields, the intra-band components dominate^, and so it 
is expected that in the limit B ^ 0, the model outlined 
here without inter-band scattering will become relevant. 

Experimental observation of an index plot with inter- 
cept consistent with 7 = in a 2DEG with small Rashba 
splitting would be a powerful confirmation of the current 
work. 



VIII. CONCLUSIONS 

We have developed a semi-classical theory of quantum 
oscillations in general particle-hole symmetry broken sys- 



tems. In particular, we have applied the formalism to 
mixed normal fermion ~ Dirac fermion systems, which 
is particularly relevant to topological insulator surface 
states. 

By properly treating the pseudo-spin magnetic mo- 
ment, we have shown that the phase offset observed in 
quantum oscillation experiments is, in general, not quan- 
tized. We found that the quantized result can be ex- 
pected in gapless systems, or in particle- hole symmetric 
systems, in agreement with previous studies. However, 
we have shown that if particle-hole symmetry is broken, 
and there is a gap, then the phase offset is not quantized. 

We have developed a protocol which allows one to de- 
termine the material properties of topological insulator 
surface states, in particular the normal fermion mass, the 
Dirac fcrmi velocity, the g-factor, and the intrinsic gap, 
by using quantum oscillation experiments. In particular, 
the effective mass cannot be naively applied to give the 
Fermi velocity, but is corrected due to Zeeman splitting 
and the energy scale mvp. Furthermore, the observed 
nonlinear index plots can be fit to a simple function at 
small fields. From this, one can read off the phase offset 
at zero field, and directly determine the intrinsic gap, or 
else regain the 'topological' result of zero phase offset, as 
expected where there is no mixing of normal and Dirac 
fermions. 

The unique interplay of the normal fermion mass and 
the Dirac Fermi velocity to create an energy scale mvp, 
together with the Zeeman or intrinsic gap, add rich sub- 
tleties to quantum oscillation experiments that allow this 
already powerful tool to probe material properties in new 
ways. 

In future studies we will focus on topologically sig- 
nificant gap phenomena, such as the topological exciton 
condensate^ and how quantum oscillation experiments 
can be used to identify and probe these topologically non- 
trivial regimes. 



Acknowledgments 

We thank ,1. Kokalj for critical readings of the 
manuscript and helpful discussions. We also thank O. 
Sushkov, A. Hamilton, and T. Li for constructive feed- 
back on the results, and especially J.N. Fuchs for provid- 
ing critical feedback of the preprint. Financial support 
was received from a UQ Postdoctoral Fellowship (ARW) 
and an Australian Research Council Discovery Project 
(RHM). 



* Electronic address: 'a. wright 7 @uq.edu. au| (2005). 

1 C.L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 226801 ^ M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 



12 



(2010) . 

^ X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011) . 

* M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, 
L.W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 
766 (2007). 

® A. Roth, C. Briine, H. Buhmann, L.W. Molenkamp, J. 
Maciejko, X.-L. Qi, and S.-C. Zhang, Science 325, 294 
(2009). 

Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. 
Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. 

C. Zhang, L R. Fisher, Z. Hussain and Z.-X. Shen, Science 
325, 178 (2009). 

D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, J. 
Osterwalder, L. Patthey, J. G. Checkelsky, N. P. Ong, A. 
V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. 
Cava and M. Z. Hasan, Nature 460, 1101 (2009). 

* D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, 
J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. 
S. Hor, R. J. Cava and M. Z. Hasan, Science 323, 5916 
(2009). 

^ T. Champel and V. P. Mineev, Philos. Mag. B 81, 55 
(2001). 

^° D.-X. Qu, Y. S. Hor, J. Xiong, R. J. Cava, and N. P. Ong, 

Science 329, 821 (2010). 
" K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 

I. Katsnelson, I. V. Grigorieva, S. V. Dubonos and A. A. 
Firsov, Nature 438, 197 (2005). 

Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

A.A. Taskin, K. Segawa, and Y. Ando, Phys. Rev. B 82, 
121302(R) (2010). 

J.G. Analytis, R.D. McDonald, S.C. Riggs, J.-H. Chu, G.S. 
Boebinger, and I.R. Fisher, Nat. Phys. 6, 960 (2010). 
Z. Ren, A.A. Taskin, S. Sasaki, K. Segawa and Y. Ando , 
Phys. Rev. B 82, 241306(R) (2010) . 

A.A. Taskin, Z. Ren, S. Sasaki, K. Segawa, and Y. Ando, 
Phys. Rev. Lett. 107, 016801 (2011). 

C. Briine, C. X. Liu, E. G. Novik, E. M. Hankiewicz, H. 
Buhmann, Y. L. Chen, X. L. Qi, Z. X. Shen, S. C. Zhang, 
and L. W. Molenkamp, Phys. Rev. Lett. 106, 126803 
(2011) . 

i« M. Veldhorst, M. Snelder, M. Hoek, T. Gang, V. K. 
Guduru, X. L. Wang, U. Zeitler, W. G. van der Wiel, A. 

A. Golubov, H. Hilgenkamp, and A. Brinkman, Nat. Mat. 

II, 417 (2012). 

L. He, F. Xiu, M. Teague, W. Jiang, Y. Fan, X. Kou, M. 
Lang, Y. Wang, G. Huang, N.-C. Yeh, and K.L. Wang, 
Nana Lett. 12, 1486 (2012). 
2" J. Xiong, Y. Luo, Y.H. Khoo, S. Jia, R. J. Cava, and N. 
P. Ong, Phys. Rev. B 86, 045314 (2012). 
J. Xiong, A.C. Petersen, D. Qu, Y.S. Hor, R.J. Cava, N.P. 
Ong, Physica E 44, 917 (2012). 

B. Sacepe, J.B. Oostinga, J. Li, A. Ubaldini, N.J.G. Couto, 

E. Giannini, and A.F. Morpugo, Nat. Comm. 2, 575 
(2011). 

A. A. Taskin, S. Sasaki, K. Segawa, and Y. Ando, Phys. 
Rev. Lett 109, 066803 (2012). 

Z. Ren, A. A. Taskin, S. Sasaki, K. Segawa, and Y. Ando, 
Phys. Rev. B 85, 155301 (2012). 

B. Seradjeh, J. Wu, and P. Phillips, Phys. Rev. Lett. 103, 
136803 (2009). 

A.A. Taskin and Y. Ando, Phys. Rev. B 84, 035301 (2011). 



^"^ B. Seradjeh, J. E. Moore, and M. Franz, Phys. Rev. Lett. 

103, 066402 (2009). 
2* G.Y. Cho and J.E. Moore, Phys. Rev. B 84, 165101 (2011). 

J.N. Fuchs, F. Piechon, M.O. Goerbig, and G. Montam- 

baux, Eur. Phys. J. B 77, 351 (2010). 
3° L.M. Roth, Phys. Rev. 145, 434 (1966). 
^1 M.V. Berry, Proc. R. Soc. Lond. A 392, 45 (1984). 

D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 

1959 (2010). 

3^ G.P. Mikitik and Y.V.Sharlai, Phys. Rev. Lett. 82, 2147 
(1999). 

L. Fu, Phys. Rev. Lett. 103, 266801 (2009). 
^5 C.-X. Liu, X.-L. Qi, H. J. Zhang, X. Dai, Z. Fang, and 

S.-C. Zhang, Phys. Rev. B 82, 045122 (2010). 
'■'^ LA. Luk'yanchuk and Y. Kopelevich, Phys. Rev. Lett. 93, 

166402 (2004). 

I. Zutic , H. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

A.P. Schnyder, S. Ryu, A. Furusaki, and A.W.W. Ludwig, 
Phys. Rev. B 78, 195125 (2008). 
»9 Z. Wang, Z.-G. Fu, S.-X. Wang, and P. Zhang, Phys. Rev. 
B 82, 085429 (2010). 

G.W. SemenofT, V. Semenoff, and F. Zhou, Phys. Rev. Lett. 
101, 087204 (2008). 

J. hinder, T. Yokoyama, and A. Sudb0, Phys. Rev. B 80, 
205401 (2009). 

Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. 
Bansil, D. Grauer, Y. S. Hor, R. J. Cava and M. Z. Hasan, 
Nat. Phys. 5, 398 (2009). 

Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, 
X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 
Zhang, I. R. Fisher, Z. Hussain, and Z.-X. Shen, Science 
325, 178 (2009). 

S. Adam, E.H. Hwang, and S. Das Sarma, Phys. Rev. B 
85, 235413 (2012). 

D. Kim, S. Cho, N.P. Butch, P. Syers, K. Kirshenbaum, S. 
Adam, J. Pagiione, M.S. Fuhrer, Nat. Phys. 8, 460 (2012). 
■"5 G.P. Mikitik and Y.V.Sharlai, Phys. Rev. B 85, 033301 
(2012). 

''■^ C.-X. Liu, H.J. Zhang, B.Yan, X.-L. Qi, T. Frauenheim, X. 
Dai, Z. Fang, and S.-C. Zhang, Phys. Rev. B 81, 041307(R) 
(2010). 

Y. Zhang, K. He, C.-Z. Chang, C.-L. Song, L.-L. Wang, X. 
Chen, J.-F. Jia, Z. Fang, X. Dai, W.-Y, Shan, S.-Q. Shen, 
Q. Niu, X.-L. Qi, S.-C. Zhang, X.-C. Ma and Q.-K. Xue, 
Nat. Phys. 6, 584 (2010). 
'^^ H.-Z. Lu, W.-Y. Shan, W. Yao, Q. Niu, and S.-Q. Shen, 
Phys. Rev. B 81, 115407 (2010). 

T. Sato, K. Segawa, K. Kosaka, S. Souma, K. Nakayama, 
K. Eto, T. Minami, Y. Ando and T. Takahashi, Nat. Phys. 
7, 840 (2011). 

X.-L. Qi, T.L. Hughes and S.-C. Zhang, Phys. Rev. B 78, 
195424 (2008). 

D. Schoenberg, Magnetic Oscillations in Metals, Cam- 
bridge University Press, Cambridge, 1984. 

'^^ J.M. Luttinger, Phys. Rev. 119, 1153 (1960). 

'^'^ P.T. Coleridge, Semicond. Sci. Technol. 5, 961 (1990). 
M. E. Raikh and T. V. Shahbazyan, Phys. Rev. B 49, 5531 
(1994). 

T.H. Sander, S.N. Holmes, J.J. Harris, D.K. Maude, and 
J.C. Portal, Phys. Rev. B 58, 13856 (1998). 



